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Abstract. All stars are born in molecular clouds, and most in giant molecular clouds (GMCs), 
which thus set the star formation activity of galaxies. We first review their observed properties, 
£SJ ■ including measures of mass surface density, E, and thus mass, M. We discuss cloud dynamics, 
' concluding most GMCs are gravitationally bound. Star formation is highly clustered within 
, GMCs, but overall is very inefficient. We compare properties of star-forming clumps with those 
/"S ■ of young stellar clusters (YSCs). The high central densities of YSCs may result via dynamical 
evolution of already-formed stars during and after star cluster formation. We discuss theoretical 
models of GMC evolution, especially addressing how turbulence is maintained, and emphasizing 
the importance of GMC collisions. We describe how feedback limits total star formation effi- 
ciency, e, in clumps. A turbulent and clumpy medium allows higher e, permitting formation of 
bound clusters even when escape speeds are less than the ionized gas sound speed. 
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' . 1. Observed Properties of Molecular Clouds 

o . 

Stars form from molecular gas, which, in the Milky Way, is mostly in GMCs, denned 
to have M > 10 4 M Q (Williams et al. 2000). Within GMCs are overdense clumps that 
may later form star clusters and cores that may form single stars or binaries. GMCs are 
normally identified by 12 CO(1-0) emission (e.g. Solomon et al. 1987 [S87]) as discrete 
, features in I — b — v space. They have a small volume fraction of the Galactic disk, but 
a significant mass fraction, ~ 1/3, of gas inside the solar circle (Wolfire et al. 2003). 

CO is dissociated by FUV photons, so CO-defmed GMCs have a protective layer of 
molecular "dark gas" (Grenier et al. 2005). For the local FUV field, dust in two layers 
on each side of the cloud provide Ay ~ 1.4 mag (van Dishoeck & Black 1988), i.e. 
S = 6.2 x lO^gcm" 2 = 3OM pc~ 2 or V H = 2.7 x 10 21 cm" 2 (we adopt n Uc /n u = 0.1, 
so mass per H, = 2.34 x 10~ 24 g). Figure 1, a diagram of £ vs. M for molecular 
clouds and star clusters, indicates this level. GMCs have larger mean £ ~ 0.02g cm~ 2 (or 
~ 100A/ Q pc -2 ) ( jjl.3p . Nevertheless, significant fractions, ~ 0.3, of total GMC molecular 
masses are predicted to be in the dark gas phase (Wolfire et al. 2010). The porous nature 
of GMCs increases FUV penetration and a large fraction, perhaps as high as 90%, of 
molecular gas can be considered to be part of Photodissociation Regions, where the FUV 
radiation field influences ionization, heating and chemistry (Hollenbach & Tielens 1999). 
The predicted equilibrium temperatures in the high extinction clumps within GMCs are 
- 10 K, corresponding to sound speeds of c s = (kT/fi) 1 / 2 = 0.188(T/10 K) 1 / 2 kms" 1 . 

1.1. Methods to Measure Mass Surface Density, S, of Molecular Clouds 

(1) From 12 CO via the "X-factor": X = iV Ha //co(i-o)- Dame et al. (2001) estimate 
X ~ (1.8 ± 0.3) x 10 20 cm -2 / (K km/s) by first finding the ratio of FIR dust emission to 
total (Hi-derived) column density, /loOjum/^H, around local, high latitude GMCs. This is 
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used together with /loo^m observed at the GMC, to predict the total iVn at its location. 
The observed contribution from atomic gas is subtracted and the remainder compared to 
^co(i-o)- Similar results are found using observed 7-ray emission (produced by cosmic 
rays interactions) to estimate total towards GMCs (e.g. Strong & Mattox 1996; 
Ackermann et al. 2011), which are then compared to HI and CO observations of these 
locations. Note, these methods allow detection of dark gas. Leroy et al. (2011) used the 
FIR emission method to measure X in 5 nearby galaxies, finding it increased by a factor 
of ~ 10 in the lowest metallicity systems, NGC 6822 & SMC — interpreted as being due 
to the increasing dominance of dark gas in these relatively dust-free environments. 

(2) From 13 CO(1-0), typically assuming its excitation temperature, T ox , is the same as 
that derived from 12 CO(1-0) (assumed to be optically thick). If the 13 CO is subthermally 
excited, then E will be underestimated. Note, Pineda et al. (2008) find that 60% of even 
the 12 CO(1-0) emission in the Perseus GMC is subthermally excited. Given an adopted 
T cx , Ni3qq can be derived, including corrections for its optical depth. Then the total 
iVn 2 follows, given an assumed abundance of 13 CO with respect to H2: e.g. Roman- 
Duval et al. (2010, RD10) assume ni2 CO /ni3 CO = 45 (Langer & Penzias 1990) and 
7ii2 CO /riH 2 = 8x 10~ 5 (Blake et al. 1987), so ni3 CO /"-H 2 = 1-78 x 10~ 6 . ni2 CO /ni3 CO is 
expected to vary with Galactic radius from ~ 40 at r = 3 kpc to ~ 70 at 8 kpc (Milam 
et al. 2005). m2 C o/?iH 2 is uncertain and may also show factor of ~ 2 variations: e.g. 
Lacy et al. (1994) find ni2 CO /n H2 ~ (3 ± 1) x 10~ 4 . Hernandez & Tan (2011) adopted 
ni2 C o/?T-i3co = 54 and ni2 CO /nn 2 = 2 x 10~ 4 so ni3 CO /nn 2 — 3.70 x 10~ 6 in their study 
of two infrared dark clouds (IRDCs). With these choices and our preferred value of /xh, 
the values of £ and M of RD10 should be scaled by a factor 0.49 (<JO]and Fig. 1). 

(3) From other molecular line tracers, e.g. C 18 0, HCO + , CS. The rarer isotopologues 
of CO, being harder to detect, tend to be observed towards higher £ clumps and cores 
within GMCs. A problem in the cold, dense conditions of IRDCs is CO freeze-out onto 
dust grains and thus its depletion from the gas phase. Widespread (extending over several 
pc) line-of-sight depletion factors, fjj, (i.e. expected gas phase abundance relative to that 
observed) of ~ 3 were seen via C 18 compared to £ derived via MIR extinction mapping 
towards an IRDC by Hernandez et al. (2011). Larger values of fjj ~ 30 were reported by 
Fontani et al. (2012) in more localized IRDC clumps via C 18 0(3-2) compared to E derived 
via FIR/submm emission. HCO + (l — 0) is often bright, but has an uncertain and likely 
variable abundance with respect to H2. It has a critical density n cr i t = 1.7x 10 5 cm -3 , but 
line trapping can lower the effective n cr it to ~ fewxl0 3 cm -3 , similar to that of optically 
thin C 18 0(1 — 0) (Evans 1999). CS can also be used to trace dense gas, especially utilizing 
the multiple higher J transitions observable from the ground (e.g. Plume et al. 1997). 

(4) From FIR/mm emission, requiring knowledge of dust temperature(s), Td, and opac- 
ity per unit gas mass, k v . Assuming optically thin transfer and black body dust emission, 
E = 4.35xl0" 3 ([5^/O]/[MJy/sr])(^/[0.01cmVg])-U?. 2 [exp(0.799T d ; i 1 5 A 1 ;2)-l]gcm- 2 , 
where A12 = A/1.2 mm and T^ig = 2^/15 K. A common choice of k v is that due to the 
moderately coagulated thin ice mantle dust model of Osscnkopf & Hcnning (1994, OH94), 
which has an opacity per unit dust mass of Ki.2mm,d = 1.056cm 2 g _1 . A gas-to-refractory- 
component-dust-mass ratio of 141 is estimated by Draine (2011), while 156 was assumed 
in the OH94 coagulation model, so «i.2mm = 6.77 x 10~ 3 cm 2 g _1 . This method probes 
all mass along the line of sight, and thus is also sensitive to dark molecular and atomic 
gas. In the Galactic plane, the method gains utility for high E clumps and cores within 
GMCs, which dominate the total line-of-sight column density. However, accounting for 
the surrounding diffuse emission is still needed (Battersby et al. 2011). 

(5) From NIR extinction mapping, using background stars (e.g. Lada et al. 2007; 
Goodman et al. 2009), probing total line-of-sight E (up to Ay ~ 25 mag). This method 
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depends on dust opacity in the NIR and the gas-to-dust mass ratio, but not T<j. It has first, 
and most accurately, been applied to nearby GMCs that sit at relatively high Galactic 
latitudes. Kainulainen et al. (2011) used this method to study clouds in the Galactic 
plane, relying on statistical binning of sources in resolution elements of ~ 30" x 30". 

(6) From MIR extinction mapping, using diffuse Galactic background emission from hot 
dust. Spitzer 8 /im images resolve to 2" and probe high-E (~ 0.5 g cm" 2 ; Ay <~ 100 mag) 
IRDCs (e.g. Butler & Tan 2009 [BT09]; 2012 [BT12]; Pcrctto & Fuller 2009). The method 
depends on K^m (BT09 use 7.5cm 2 g" 1 [OH94]), but not T<j. Allowance is made for fore- 
ground emission, best measured by finding "saturated" intensities towards independent, 
optically thick cores (BT12). Only differences in E relative to local surroundings are 
probed, so the method is insensitive to low-E environs of IRDCs. This limitation is ad- 
dressed by combining NIR & MIR extinction maps (Kainulainen & Tan 2012 [KT12]). 

Once E is measured, cloud mass, M, is found from projected angular area and an 
estimate of distance, typically kinematic for Galactic clouds. In the inner Galaxy, the 
near/far degeneracy may be broken by looking for HI self-absorption. For IRDCs, the 
near kinematic distance is often assumed. Kinematic distances are uncertain due to non- 
circular cloud motions caused by their ^5—10 km s" 1 disk-plane velocity dispersion or 
streaming due to spiral arms or bars. Parallax distances to masers associated with star 
formation are available for some sources (see, e.g. Reid et al. 2009; Foster et al. 2012). 

1.2. Virial Equilibrium of Molecular Clouds 

The ID internal velocity dispersion, a, of clouds is measured from their line widths. 
This allows assessment of cloud kinetic energy, T = (3/2)Mct 2 , and comparison to 
gravitational energy, W — — (3/5)aGM 2 / R, where a is a factor, <~ 0(1), that ac- 
counts for nonspherical and nonuniform density distributions (Bcrtoldi & McKcc 1992 
[BM92]). BM92 find the effects of nonsphericity are quite small, so we neglect them 
here. For a power-law density distribution p oc r~ kp , a = (1 — fc p /3)/(l — 2fc p /5) = 
10/9, 5/4, 5/3 for k p — 1, 3/2, 2, respectively. BM92 define the virial parameter as a v ir = 
5a 2 R/(GM) = 2aT/\W\. A cloud in virial equilibrium satisfies 2(T - To) + M + W = 0, 
where To = (3/2)P ^ (Po IS the surface pressure around the cloud volume, V) and 
M = [l/(87r)] J (B 2 — B 2 )dV (B is the magnetic field far from the cloud and the inte- 
gral extends beyond the cloud surface over regions where the B-ficld suffers distortions) . 

Molecular clouds are magnetized (Crutcher 2012), with observed values being consis- 
tent with a uniform distribution of total field strengths from to a maximum, -Btot.o = 
IQfiG for n H < 300cm" 3 , and B TO t,o = 10(n H /300 cm" 3 ) 065 for n H > 300 cm" 3 . How- 
ever, in general we do not have accurate estimates of total field strengths in and around 
a given CMC or clump, making it difficult to assess its contribution to virial balance. 

For a cloud satisfying the equilibrium virial equation, a v i r = a(l + (27o — .M)/|W|) = 
a(l — Po/-P) _1 (l — A1/IT4 7 !), where P = (3/2)pcr 2 is the volume-averaged cloud pressure. 
For negligible surface pressure (P a — > 0) and magnetic support (M — > 0), a v ; r = a, i.e. 
T = \ W\/2. Then, a gravitationally bound cloud has a va < 2. If a cloud has P = P/2, 
a vil is raised by a factor of 2. If M = \W\/2, lowered by a factor of 2. 

1.3. Giant Molecular Clouds 
For 12 CO-dcfincd_GMCs, S87 measured E - 170 M Q pc" 2 . However, Heyer et al. (2009, 
[H09]) measured E ~ 42 M Q pc" 2 for these same structures using 13 CO. As discussed by 
H09, the true E is likely to be lower than the S87 estimate due to incomplete sampling 
of the original 12 CO maps, but the 13 CO value may be underestimated because of sub- 
thermal excitation. We indicate the properties of Galactic 12 CO-defined GMCs in Fig. 1 
with E - 100 M© pc" 2 and a mass range from 10 4 M Q to - 6 x 10 6 M© (Williams & 
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McKee 1997). Several outer-Galaxy GMCs and their clumps were studied by Heyer et 
al. (2001), who found they had E ~ 10 M pc" 2 . Bolatto et al. (2008) studied GMCs in 
nearby galaxies, finding they generally had properties similar to Galactic GMCs. How- 
ever, GMCs may differ in more extreme environments: e.g. 14 GMCs observed in the 
dwarf starburst Henize 2-10 have E ~ 400 M Q pc" 2 (Santangelo et al. 2009 [S09]). 

The 13 CO-defined Galactic GMCs and clumps of RD10 (with E & M scaled by 0.49, 
[EH) are also shown in Fig. 1. For clouds with M > 10 4 M Q , E = 124M©pc~ 2 , including 
dark gas. For this sample, considering a V ir reported by RD10 (increased by factors 2.02 
to account for the mass scaling [ §1.1| and 1.28 to match the BM92 definition, for total 
of 2.60; note, we neglect contribution of dark gas to a V ir) we find a V ir = 1-09. Also, 89% 
of these 249 GMCs (with 95% of the mass) have a V ir < 2, i.e. are gravitationally bound. 
Similar results are obtained for 12 CO-defined GMCs, if one allows for a factor ~2 increase 
in H09 masses (above) and for the extragalactic GMCs studied by Bolatto et al. (2008). 

GMCs that are gravitationally bound and near virial equilibrium explain H09's result of 
a modified line- width vs. size relation: a oc (Ei?)°- 5 , which follows if M ~ M v - n ~ 5a 2 R/G 
and E = M/(ttR 2 ). Note, this may also explain the higher normalization of the relation 
seen in high-E clumps (e.g. Plume et al. 1997), if they are also self-gravitating. 

Thus we conclude that GMCs are gravitationally bound and have configurations that 
are consistent with being in approximate gravitational virial equilibrium (a V i r ~ 1). 

1.4. Clumps and Cores 

Barnes et al. (2011 [Bll]) derived the properties of a complete census of ^300 HCO + 
clumps in a 20° x 6° region of the Galactic plane (Fig. 1). While the absolute properties of 
individual sources are uncertain due to the assumed HCO + abundance and T cx , relative 
properties should be more accurate. The most extreme (in M & E) clump in this sample 
is G286 (Barnes et al. 2010 [B10]), which also contains an embedded, forming cluster. 

Ginsburg et al. (2012 [G12]) presented properties of a census of massive (> 10 4 Mq) 
clumps found by 1.1 mm emission in the Bolocam Galactic Plane Survey (BGPS). In 
Fig. 1 we have scaled their (free-free-subtracted, Td = 20 K) masses by a factor of 1.47 
to account for our preferred choice of fci.i mro = 7.78 x 10 3 cm 2 g _1 (OH94). Longmore et 
al. (2012 [L12]) reported a comparable, cold, 1.3 x 1O 5 M , 2.8 pc-radius Galactic center 
clump, dubbed the "Brick" . It shows little evidence of star formation activity. 

Cold, high-E regions of GMCs reveal themselves as IRDCs and 10 such clouds are 
shown in Fig. 1, with masses derived from the combined NIR+MIR extinction mapping 
method of KT12. We also show 42 core/clumps selected from these IRDCs by BT12, 
(properties derived for best-fitting singular polytropic sphere model). The structure is 
resolved interior to these radii down to 2" (~ 0.04 pc at 4 kpc). The mean value of k p = 
1.1 ± 0.25. If account is taken of the surrounding clump envelope, then k p = 1.6 ± 0.28. 
Some of these objects, as well as others reported by Bontemps et al. (2010), appear to be 
consistent in their properties with massive starless cores invoked in the Turbulent Core 
Model of massive star formation (MSF) (McKee & Tan 2003 [MT03]). 

Mueller et al. (2002 [M02]) studied dense cores and clumps, selected by water maser 
activity to harbor MSF. For 31 sources, they derived the density and temperature struc- 
ture by ID radiative transfer modeling of the IR/submm emission, finding k p ~ 1.8 ±0.4 
(however, these results may be affected by the presence of bipolar outflow cavities, Zhang 
& Tan 2011). For consistency, in Fig. 1 we have multiplied the M02 values of E and M 
by a factor of 1.56, and we show their properties at a scale where «h~3x 10 4 cm~ 3 . 

Most lower M and E 13 CO and HCO + clumps appear pressure confined (e.g. BM92, 
Bll), i.e. a v i r S> 1 and their internal pressure is influenced more by the ambient pressure 
rather than by self-gravitating weight. Note, the clumps may still satisfy the equilibrium 
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virial equation. Higher-E clumps appear to be nearer gravitational virial equilibrium (i.e. 
internal pressure set by self- gravity to be ~ G£ 2 ). This has been found in IRDC filaments 
(Hernandez et al. 2012), clumps (KT12) and cores (Tan et al., in prep.). In G286, BIO 
report a mass infall rate 3.4 x 10 -2 M Q yr" 1 , highest in the Bll sample, but still ~10 
times less than the rate expected under global free-fall collapse. 



o 
o 



B 

o 
w 



o ; 



o 
d 



\j \ i i 1 1 1 1 1 1 — i i.-i 1 1 iij 1 — i..n 1 1 1 1 1 1 — i i.-i 1 1 1 1 1 i^\,.n i in 1 — i i-i 1 1 1 1 



% 



or- 



Turbulent Core M-odel 
(M=60M o , E cl =l g cpi'-') ..- 

L- IRDC * .../•"'" 
Core/Glumps ro t e ' 



rches.. 



+ + 



HC0 + 



M82 



v.* • ^ 



... -;J<GC5g53 SSC 



SSCs- 



JjgC1569-A^ 

fir- 




y^-.Charn-ps/ V 

' T 

ouds 1 



J I I t'l llll I -I I I.Hiil | I I M I III 



_ CO^ clouds J . i^i - 

gas ' \ 

_l i i i-i-riil I i i t'l 1 1 1 1 i - 1 i i.nii 



10 



100 



1000 



10 4 



10 5 



10 6 



10 7 



M (M ) 

Figure 1. E-M Diagram of Molecular Clouds & Young Star Clusters. Mass surface density, 
E = M/(nR 2 ), is plotted vs. mass, M. Dotted lines of constant radius, R, H number density, nn 
(or free-fall time, t B = [37r/(32Gp)] 1/2 ), and escape speed, v csc = (10/« v ir) c, are shown. The 
minimum E for local CO-emitting clouds due to FUV radiation is ~ 30 Mq pc -2 (dashed line). 
Below this molecular "dark gas" can exist. The following objects are plotted (see text for more 
details). Typical 12 CO-defined GMCs have E ~ 100 Mq pc -2 , although denser examples have 
been found in Henize 2-10 (S09). The 13 CO-defined clouds of RD10 are indicated, along with 
HCO+ clumps of Bll, including G286 (B10). 10 IRDCs (KT12) and their internal core/clumps 
(BT12) are shown, the latter overlapping with Massive Star-Forming core/clumps [M02]. Along 
with G286, the BGPS clumps (G12) and the Galactic Center Brick (L12) are some of the most 
massive, high-E gas clumps known. Clumps may give rise to young star clusters (see text for 
refs.), like the ONC and NGC 3603 (radial structure shown from core to half-mass, R1/2, to outer 
radius), or even more massive examples, e.g. Westerlund 1, Arches, Quintuplet (shown at R\/2)- 
Example clusters in the LMC (R136) & SMC (NGC 346) show a wide range of E. More massive, 
Super Star Clusters are found in some dwarf irregular galaxies, NGC 1569 & 5253, and starburst 
galaxy M82. Grey squares indicate the GMCs simulated by VLBT12 from beginning (large) to 
end (small) of a 10 Myr period, including merger and fragmentation of clouds. Grey circles, A, B, 
C, show TM01 feedback model clumps, yielding star formation efficiencies of e = 0.33, 0.48, 0.50, 
respectively. Grey pentagon shows initial condition of a fiducial MSF turbulent core of MT03. 
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1.5. Star Formation: Embedded and Young Star Clusters (YSCs) 

The star formation efficiency per ts of GMCs and clumps is puny, es ~ 2% (Zuckerman & 
Evans 1974; Krumholz & Tan 2007). It is unclear whether this is due mostly to support 
by turbulence (Krumholz & McKee 2005) or magnetic fields, although evidence for a 
threshold E ~ 120 M pc -2 (Ay ~ 8 mag) for star formation (Lada et al. 2010) may 
indicate the importance of the latter (McKee 1989). Star formation is concentrated in 
clumps in GMCs, where total efficiencies e > 0.3 can be achieved, as evidenced by bound 
YSCs. A ratio egf/e <C 1 implies cluster formation times, t{ olm ^> tg (Tan et al. 2006). 

Embedded star clusters provide a snapshot of the star cluster formation process. YSCs 
may also retain memory of the physical conditions of the gas clumps from which they 
formed, although modified by mass segregation, cluster core collapse, expansion due to 
gas massloss, and tidal stripping (see Portegies-Zwart et al. 2010 [PZMG10]). In Fig. 1 
we show the Orion Nebula Cluster (ONC), with density structure from Hillenbrand & 
Hartmann (1998) with M and E estimated at the core, half-mass, R1/2, and total radii. 
These masses are virial mass estimates, and so receive contribution from both stars and 
gas. Similarly, we show NGC 3603 — a significantly more massive cluster (Harayama et 
al. 2008). The most massive known Galactic YSC is Westerlund 1, with dynamical mass 
- 6.3±|:f x 10 4 M Q and R 1/2 =0.86 pc (Mengel & Tacconi-Garman 2007). Harfst et al. 
(2010) modeled the expected initial conditions of the Arches cluster, near the Galactic 
Center, to match its observed stellar density distribution 2.5 Myr later. Their preferred 
best-fit model has an initial mass of (4.9 ± 0.8) x 10 4 M Q (IMF uncertainties introduce 
larger errors) and a virial radius of 0.70 pc (i.e. R1/2 — 0.44 pc). We show this half-mass 
condition in Fig. 1. The Quintuplet cluster currently has about 6,000 Mq inside a radius 
of 0.5 pc (Hufimann et al. 2012), but being ^4 Myr old, it is more difficult to reconstruct 
its initial conditions: we adopt M = 4x 10 4 M Q , R 1/2 =2 pc (Figer et al. 2006; PZMG10). 

In the 30 Doradus region of the LMC, R136 has M ~ 6 x 10 4 M (Andersen et al. 2009) 
and R1/2 =1.7 pc (Hunter et al. 1995), similar to the dynamical mass estimate of Henault- 
Brunet et al. (2012). NGC 346 is a much lower density star cluster (or association) 
in the SMC, with M = 3.9 x 10 5 Af© and R 1/2 = 9 pc. More massive, "Super Star 
Clusters" (SSCs) have been found in, e.g., the dwarf irregular galaxies NGC 1569 (cluster 
A with M = 1.6 x 10 6 M 0) R 1/2 = 2.3 pc, Anders et al. 2004, PZMG10; and B with 
M = 4.4 x 10 4 M , Rx/2 = 2.2 pc, Larsen et al. 2008) and NGC 5253 (M ~ 5x 10 5 M Q , 
R x /2 ~ 0.6 pc, Turner & Beck 2004) and in M82 (McCrady & Graham 2007) (Fig. 1). 

A comparison of YSCs and gas clumps in the E — M diagram (Fig. 1) suggests relatively 
massive starless (or at least gas-dominated) clumps do form, albeit rarely, that may be 
progenitors of the most massive Galactic clusters. These can have fin ~ 10 5 cm -3 and 
tff ~ 10 5 yr. Star formation appears inhibited during assembly of these massive clumps. 
We suspect this is due to support by magnetic fields rather than turbulence, which 
would quickly decay (fj2]). It is not yet clear if ~ 10 6 M starless clumps can form to 
be progenitors of SSCs. Central core densities of many YSCs exceed what is seen in gas 
clumps and may be the result of dynamical evolution of the cluster both during the star 
cluster formation process (especially if tf orm 3> tg) and after star formation has ceased. 

2. Theoretical Models of Molecular Clouds 

GMCs and clumps exhibit density structures on a wide range of scales, all the way down 
to those of forming stars. They contain motions, initially imprinted from the surroundings 
or due to self-gravity, and later due to feedback processes, that are many times the sound 
speed, leading to strong shocks. They are magnetized. They have poorly-defined bound- 
aries, merging continuously into the surrounding galactic environment. These features 
make theoretical and numerical modeling of molecular clouds a challenging problem. 
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2.1. Turbulence and GMC Evolution 



GMCs velocity dispersions are 3> c s . Supersonic turbulence decays in ~ 1 dynamical 
time, idyn — R/cr (Stone ct al. 1998; Mac Low et al. 1998) (for a V i r — l f iff — 0.5idyn)- 
Maintenance of turbulence is a constraint on models of GMC formation and evolution. 
McKee & Ostriker (2007 [MO07]) reviewed two conceptual frameworks: (1) GMCs are 
dynamic, transient and largely unbound, with turbulence driven by large-scale colliding 
atomic flows that form the clouds (e.g. Heitsch et al. 2005; Vazquez-Semadeni et al. 2011). 
These models do not explain why most GMCs are bound with a v ; r ~ 1 f £|1.3[) . Fast flows 
of HI around GMCs needed to form them quickly in ~ td yn have not been observed. Nor 
is it clear if such models are globally consistent with the relatively high mass fractions of 
gas in GMCs, at least inside the solar circle. (2) GMCs form by large-scale gravitational 
instabilities and are thus gravitationally bound with a V ir ~ 1- Turbulence is maintained 
by contraction and then, later, star formation feedback (e.g. Goldbaum et al. 2011). 

Here we present a third paradigm: GMC evolution under frequent cloud collisions (Tan 
2000; Tasker & Tan 2009 [TT09]; Tan & Shaske, in prep.). In a galactic disk with Toomre 
Q ~ 1 and significant fraction, /gmC ~ 1/2, of ISM gas mass in gravitationally bound 
clouds, (centers of which appear as GMCs), clouds suffer collisions on timescale, 



icon — [26' (tt - dv c 



2 /dr)N A r 2 t f G ] 



(2.1) 



where b' = b/r t ~ 1.6 is typical impact parameter of collisions (Gammie et al. 1991 [G91]), 
= w C i rc /r = 27r/£ or b, w c irc is orbital speed at galactocentric radius r, Ma is number of 
GMCs per unit area, r t = (1 - /3)~ 2/3 (2M/M gal ) 1 / 3 r is GMC tidal radius of GMC at 
location of given /3 = dhw c i vc /dlnr and interior galactic mass M ga i, and fo ~ 0.5 is frac- 
tion of strong gravitational encounters that lead to collision. For Ma — fcMC^gas/M and 
a gas ~ (GAf K ) 1 /3(i. - 1J/3) (G91), we have f co n ~ 0.15Qi orb /(67GMc/ G [l - 0.7/3]) 
0.19/Q ] J Ic i or b, where the last evaluation is for constant u c iro i-e. (3 — 0. TT09 found 
approximate agreement with this result in their simulation of GMCs in a galactic disk. 
In Fig. 2a we evaluate eq. flUT) for each TT09 GMC from t = 225 - 275 Myr, and com- 
pare the predicted mean f co \\ = t C oii/^orb with that measured in the simulation. Good 
agreement is found for both the entire GMC sample and those with M > 10 6 Mq. 

We balance the momentum injection rate due to cloud collisions, pec, with the dis- 
sipation rate within the clouds, pdiss, to predict typical GMC density and E. We have 



Pec 



where the relative velocity v 



ar t il, for which v s ^ rt is the 




Figure 2. (a) Left: GMC-GMC collision time, t co n, relative to local orbital time, i or b vs. 
in the TT09 simulation for all (solid line) and M > 10 6 Mq (dashed line) GMCs. Solid/open 
squares show the analytic estimate (eq. I2.1|l for these respective populations (vertical lines show 
dispersion of the averages sampled 50 times at 1 Myr intervals, (b) Right: Volume density 
rendering (hh ~ 100 cm -3 - blue; 10 5 cm -3 - red) of a 0.5 pc-resolution simulation of 1 kpc 
patch of a galactic disk (VLBT12), extracted from the TT09 simulation and evolved for 10 Myr. 
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shear velocity at b — r t and a ~ 1.9 is expected for typical b = 1.6r t added in quadra- 
ture with (Tg as ~ w s ,r t - Thus, pec — (a/[7r/coii])(2G) 1 / 3 (Mw C i rc /r) 4 / 3 . For turbulent dis- 
sipation we have p d i ss = ~p/(a diss t s ) = ~V3Ma/(a diss t s ) = -2.19aV 2 GME/a d i ss . 
The turbulent energy dissipation rate has been measured in numerical simulations to be 
e d iss = (l/2)(E/E)[l /(V3a)} ~ 0.6 (see review by MO07), so that a diss ~ 3.85 (in a 
cloud with a V ir — 1 and Iq = 2R) . Thus the equilibrium value of E for GMCs is 



This fiducial value is similar to E ~ 300 M pc~ 2 of the TT09 GMCs, suggesting 
momentum injection from collisions is maintaining their turbulence and setting their 
density. Several GMCs extracted from the TT09 simulation have been evolved to higher, 
0.5 pc, resolution by Van Loo et al. (2012 [VLBT12]) (Fig. 2b). Dense clumps and 
filaments form, some reminiscent of IRDCs such as the ~ 100 pc-long "Nessie" filament 
(Jackson et al. 2010), but overall Es of the GMCs remain quite constant (Fig. 1), with 
clouds supported by a turbulent energy cascade from larger scales, including via collisions. 

Evaluating eq. (|2 .2[) for Galactic conditions, where /gmc, / C oU an d M are somewhat 
smaller than in the TT09 simulation, then we obtain E eq closer to the ~ 100 M Q pc -2 
of observed GMCs. In addition to helping to maintaining turbulence, relatively frequent 
GMC collisions produce a wide spread in rotation directions (Dobbs 2008; TT09), split 
roughly evenly between pro and retrograde with respect to galactic rotation, as is seen 
in Galactic (Koda et al. 2006; Imara k Blitz 2011) and M33 (Imara et al. 2011) GMCs. 



How do so many O stars form together in the small volume of a massive YSC in the 
face of their mutual feedback of winds, ionizing photons and radiation pressure? Tan & 
McKee (2001 [TM01]; 2004), modeled feedback on clumps of given M & E, idealized as 
a population of cores in an intercore medium. To approximate turbulent motions, cores 
orbited freely in the clump potential. Stars were formed at constant rate, eg = 0.02, with a 
near-Salpeter IMF, and placed at clump center. Wind, ionization and radiation pressure 
feedback were derived from Starburst99 models. The spherical average extent of the 
HII region was calculated balancing ionizing flux with recombinations, dust absorption 
and core shadowing. Cores exposed to ionizing flux were imploded with the magnetized 
radiation driven implosion models of Bcrtoldi & McKee (1990). Core dynamics included 
the photoevaporative rocket effect, as well as radiation pressure. Advection of neutral 
cores into the HII region and subsequent disruption maintained high ionized gas densities, 
thus confining the HII region. A momentum conserving stellar wind bubble swept-out 
the HII region center, its boundary set by balancing wind ram pressure with HII region 
thermal pressure. Star formation was stopped once the HII region neared the edge of the 
clump or after 3 Myr, at the onset of supernovae. Three example clumps were considered, 
A, B, C (Fig. 1), with M = 4 X 10 3 ,4 x 10 4 ,4 x 10 5 M and E = 1.0, 2.2, 4.6 g cm" 2 , 
resulting in e = 0.33, 0.48, 0.50, respectively. Feedback confinement and e high enough to 
form a bound cluster occurred even when v csc < lOkms^ 1 , the ionized gas sound speed. 

Fall et al. (2010) presented a simpler analytic model that considers energy and mo- 
mentum driven feedback on a clump with a smooth radial density profile. Numerical 
simulations of feedback interaction with clumpy and/or turbulent gas have been pre- 
sented for: ionization by, e.g., Henney et al. (2009), Gritschneder et al. (2010), and Dale 
et al. (2012); for wind/supernova shocks by, e.g., Pittard et al. (2010); and for radiation 
pressure by, e.g., Krumholz & Thompson (2012). However, the challenge remains of com- 
bining and fully resolving all the processes thought important for star cluster formation. 
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